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Abstract 

Background: Heterochronic shifts during ontogeny can result in adaptively important innovations and might be 
initiated by simple developmental switches. Understanding the nature of these developmental events can provide 
insights into fundamental molecular mechanisms of evolutionary change. Fishes from the Suborder Belonoidei 
display a vast array of extreme craniofacial morphologies that appear to have arisen through a series of heterochronic 
shifts. We performed a molecular heterochrony study, comparing postembryonic jaw development in representatives 
of the Suborder Belonoidei, the halfbeak Dermogenys pusillo (where the lower jaw is considerably elongated compared 
to the upper jaw) and the needlefish Belone belone (where both jaws are elongated), to a representative of their sister 
group the Suborder Adrianichthyoidei, the medaka Oryzias latipes, which has retained the ancestral morphology. 

Results: Early in development, the lower jaw displays accelerated growth both in needlefish and halfbeak compared to 
medaka, and secondary acceleration of the upper jaw is seen in needlefish later in their development, representing a 
case of mosaic heterochrony. We identified toothless extensions of the dentaries as innovations of Belonoid fishes and 
the source of heterochronic growth. The molecular basis of growth heterochronies in the Belonoidei was examined 
through comparing expression of skeletogenic genes during development of halfbeak and medaka. The calmodulin 
paralogue calm! was identified as a potential regulator of jaw length in halfbeak as its expression gradually increases in 
the lower jaw, but not the upper jaw, in a pattern that matches its outgrowth. Moreover, medaka displays equal 
expression of calm! in the upper and lower jaws, consistent with the lack of jaw outgrowth in this species. 

Conclusions: Heterochronic shifts in jaw growth have occurred repeatedly during the evolution of Belonoid fishes and 
we identify toothless extensions of the dentaries as an important innovation of this group. Our results suggest that 
calm 1 contributes to jaw heterochrony in halfbeak, potentially driving further heterochronic shifts in jaw growth across 
the Suborder Belonoidei, such as the upper jaw acceleration observed in needlefish. 

Keywords: Heterochrony, Morphological innovation. Jaw development, Beloniformes, Gavin de Beer, Dermogenys, 
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Background 

One of the major goals of evolutionary biology is to 
understand the developmental and genetic origins of 
morphological novelties that are of adaptive relevance [1]. 
Already classic work by Gavin De Beer advanced the idea 
that some novelties arise as the result of heterochronic 
shifts, which involve an alteration in the timing or rate of 
ontogenetic events during evolution, rather than having a 
completely de novo origin [2-8]. Although categorizations 
of heterochronic shifts have been devised (Additional file 1: 
Figure SI), they essentially involve - compared to the 
ancestral condition - either an increased growth rate or 
the addition of ontogenetic steps (peramorphosis), or the 
retardation of growth or a reduced number of ontogenetic 
steps (paedomorphosis) [9]. While heterochronic shifts 
in growth can result in dramatic alterations to adult 
morphology, due to the modular nature of development, 
such changes can be induced by simple developmental 
switches [10]. 

Increasingly, molecular developmental approaches are 
becoming incorporated into studies of heterochrony since 
they can shed light on the genetic bases of the heterochro- 
nic shifts and the regulatory changes that may underlie 
them [4,7,8,11]. By applying a comparative framework, 
molecular studies of heterochrony have the potential to 
identify not only the ontogenetic stage, structure and 
magnitude of growth of morphological alterations, but 
also their genetic bases. This approach involves develop- 
mental comparisons between members of different line- 
ages to their sister lineages that retained the ancestral 
condition (for examples see [8,12-15]). Heterochronic 
growth may be underlain by various developmental 
mechanisms that alter the proliferation, differentiation 
or apoptosis of various cell populations and as such, identi- 
fying its genetic basis can be aided by a fine-scaled develop- 
mental approach [4,16]. Such developmental investigations 
have been successfully utilized in studies in cichlid fish and 
Galapagos finches, which interestingly, point at similar 
molecular drivers of differential proliferation in cellu- 
lar condensations that contribute to shape differences 
in adult jaws and beaks respectively (that is, bmp4 and 
calmodulin) [17-22]. 

For more than 100 years, due to their striking jaw 
morphologies beloniform fishes have been subjected to 
ontogenetic studies and, even more importantly, inspired 
giants in the field such as Severtzov, De Beer, and Gould 
to devise the theory of heterochrony ([2,23]; [24-26] cited 
in [27]). In fact the beloniform example played an import- 
ant role in the incorporation of Darwinian thinking into 
developmental biology and the founding of the field of 
evolutionary morphology (reviewed in [3,28]). In Severtzov s 
classic treatment [29] he sorted out some of the confusions 
and misinterpretations of von Baers law and from the 
misguided interpretation of Haeckel. 



The Order Beloniformes includes the Suborder Belonoidei, 
which contains such diverse representatives as the nee- 
dlefishes (Belonidae), halfbeaks (Hemiramphidae), flying 
fishes (Exocoetidae) and sauries (Scomberosocidae). The 
single other Suborder in the Order Beloniformes is the 
Adrianichthyoidei, the ricefishes, which includes medaka, 
a developmental model system. The families in the Sub- 
order Belonoidei differ remarkably in the relative lengths 
of their upper and lower jaws. Halfbeaks, as the name 
suggests, have elongated lower jaws and relatively short, 
non-pro trusible upper jaws whereas needlefish have elon- 
gated upper and lower jaws of equal length. Interestingly, 
most needlefish pass through a 'halfbeak' stage as juve- 
niles, which led to the hypothesis that halfbeaks might be 
paedomorphic derivatives of the needlefish ([26] cited in 
[2,27]). This hypothesis was called into question however, 
by a molecular phylogeny of the Beloniformes, which 
formed the basis for a reconstruction of ontogenetic 
transitions in jaw length (Figure lA) [27,30]. The secondary 
growth of the upper jaw of the needlefish would therefore 
best be interpreted as peramorphosis (specifically, hyper- 
morphosis; see Additional file 1: Figure SI) rather than 
paedomorphosis. 

Within the Belonoidei the 'halfbeak' morphology appears 
to represent the ancestral state, from which all other jaw 
morphologies are derived (including the secondary resorp- 
tion of an initially elongated lower jaw in flying fish, as 
would be implied by the phylogeny). Interestingly, the 
members of the Suborder Adrianichthyoidei never display 
the 'halfbeak' morphology and predominantly lack elong- 
ate jaws (with the exception of some members of the 
genus Adrianichthys, which display jaw elongation that is 
either isometric {A, poptae) or even longer in the upper 
jaw {A, krutyi and A roseni)) [31]. The flexibility of their ra- 
tio in jaw length appears to have promoted the evolutionary 
success of the Belonoidei compared to the Adrianichthyoi- 
dei, which include 240 and 33 species respectively, based 
on current estimates from FishBase (http://www.fishbase. 
org). This suggests an enhanced degree of evolvability, 
defined by Hansen as 'the ability of the genetic system to 
produce and maintain potentially adaptive genetic variants' 
[32]. Pronounced variation of jaw ratios in the Belonoidei, 
driven by growth heterochronies, appears to have promoted 
the exploitation of novel ecological niches over brief evolu- 
tionary timescales, which might explain their considerably 
greater species richness in comparison to their sister group. 

While previous research has provided strong evidence 
that heterochrony has played a significant role in the 
evolution of beloniform fishes, it has not addressed the 
developmental or genetic bases of it. Importantly, previous 
studies on beloniform fishes did not compare their jaw de- 
velopment to that of medaka, a sister group [33] that has 
retained the ancestral condition of short upper and lower 
jaws, the character state that is observed in most fishes. 
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Figure 1 Heterochrony has contributed significantly to craniofacial diversity in the Belonoidei. (A) Simplified phylogeny based on [27] 
showing the relationships between representative taxa of the Belonoidei including needlefish and halfbeak, and medaka, which belong to their 
sister group the Adrianichthyoidei. (B, C) Developmental heterochrony was assessed in the jaws of these three species by plotting shape, represented by 
the ratio of lower jaw (Ij) to upper jaw (uj) length (lj:uj), against relative age. The term 'onset age' refers to the commencement of skeletal ossification, while 
'offset age' refers to the completion of ossification. Medaka retained the ancestral condition, maintaining an lj:uj of 1.0. Conversely, needlefish and halfbeak 
display acceleration of Ij growth relative to the uj (indicated by lj:uj > 1.0). This ratio is maintained to adulthood in halfbeak. Additionally, needlefish 
display a secondary acceleration of uj growth, resulting in a ratio of approximately 1.0, which is then maintained until adulthood. (D, E) In addition to 
the aforementioned heterochronic shifts, the uj and Ij also display allometry in their growth, as indicated by plotting their ontogenetic trajectories. 
(C E) 95% confidence intervals are indicated for both the heterochrony plot and ontogenetic trajectory respectively. 



Here we present a molecular developmental study of jaw 
elongation in the Beloniformes, with a comparative ana- 
lysis of candidate gene expression during jaw outgrowth 
in halfbeak and medaka. The candidate genes include two 
markers of cartilage and bone development, sox9 and 
runx2 respectively, and genes associated with differential 
cellular proliferation in skeletal primordia, including bmp4 
and calmodulin, and hmp2. This approach enables a more 
informed interpretation of the significance of any changes 
in gene expression within a developmental framework. 

Methods 

Collection of juveniles 

Adults of Dermogenys pusilla, a live-bearing species, were 
obtained through the aquarium trade in 2008, raised at a 
low density in well-planted tanks and fed flake food and 



frozen Drosophila melanogaster. Pregnant females were 
isolated in 40-litre tanks to give birth, to avoid canniba- 
lization of the clutch by conspecifics. Larvae were collected 
at the required ages and euthanized with MS222, accord- 
ing to animal ethics requirements (Regierungsprasidiums 
Freiburg, Germany). Specimens used for skeletal analysis 
were fixed with 4% paraformaldehyde in 0.1 M PBS at 
4°C, and stepped into 100% MeOH. Larval qRT-PCR 
samples were collected in the same way and stored in 
RNAlater (Qiagen, Stockach, Germany) until further pro- 
cessing. For the pre-birth stage, pregnant females were eu- 
thanized in MS222 and embryos were dissected from the 
uterus and stored in RNA later. Embryos of Oryzias latipes 
were collected using methods described in [34] and raised 
at 25°C with a few drops of methylene blue. Samples were 
fixed using the same methods as described for D, pusilla. 
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Histochemistry and measurements 

Skeletal staining was performed on fixed specimens using 
a modified alcian blue and alizarin red method, which uses 
acid-free alcian blue containing 60 mM MgCl2 [35] . Three 
specimens for each time point were photographed with a 
Zeiss Stemi SV 11 microscope (Zeiss, Munich, Germany). 
Brightness and contrast were adjusted with Photoshop 
CS4 (Adobe Creative Cloud) and measurements of key 
features were made using ImageJ. Measurements of upper 
and lower jaw lengths were taken from the laterally 
oriented specimens, calculating the distance from the 
point of attachment between the maxilla and dentary, to 
the anterior tips of the premaxilla and dentary (toothless 
extensions in the case of the halfbeak and needlefish). We 
then analysed both heterochrony and allometry, using 
methods modified from Alberch et al [9] and Boughton 
et al, [23]. The shape variable selected for our hetero- 
chrony analysis was the ratio between lower and upper 
jaws, as jaw ratio is considered to be an important shape 
measure, as this measure appears to relate to trophic 
niche according to previous studies [23,27,30]. We scaled 
the ages of all three species, similar to Alberch et al [9] 
and Gould [36] such that onset age was birth or hatching 
for halfbeak and needlefish respectively and medaka were 
age-matched on the basis of skeletal ossification. Offset 
age was at completion of skeletal ossification for all three 
species. Allometry analyses (ontogenetic scaling) plotted 
growth of uj against Ij. Additionally, 10 (iM sections were 
cut for paraffin-embedded specimens and stained using 
acid-free alcian blue and alizarin red. 

Isolation of genes of interest 

Candidate genes were isolated from cDNA of a single D. 
pusilla larva shortly after birth. Fragments of candidate 
genes were amplified using primers designed to alignments 
of teleost sequences obtained from GenBank (Additional 
file 2: Table SI) and cloned into either PGEM-T (Promega, 
Mannheim, Germany) or PCR-4 (Invitrogen, Darmstadt, 
Germany) vectors. RACE PCR and cloning was used to 
obtain 3' UTR sequences. To determine the orthology 
of our candidate genes, we constructed maximum likeli- 
hood (ML) trees. Sequences were aligned to other verte- 
brate orthologues and paralogues in Jalview [37] using 
the MAFFT algorithm. After removal of gaps, ML trees 
were generated in PhyML [38], with 1,000 bootstrap 
replicates. Calmodulin orthologues were obtained from 
medaka through the GenBank database. To determine 
whether the genes of interest had undergone positive 
selection, we calculated the ratio between the number of 
non-synonymous substitutions (dN) per non-synonymous 
site versus the number of synonymous substitutions 
per synonymous site (dS) [39]. To achieve this we 
compared the halfbeak coding sequences to other 
teleost representatives, using ynOO in the PAML toolkit 



(v. 4.4) [40], selecting the Yang and Nielsen substitution 
model [41]. 

qRT-PCR of candidate genes in the jaws of Dermogenys 

The upper and lower jaws of specimens were excised under 
a dissecting microscope at the level of the nares (roughly 
corresponding to the point at which the maxilla is attached 
to the dentary). Upper and lower jaw samples include 
skeletal elements listed in Additional file 3: Table S2. RNA 
was extracted with the RNeasy Minikit (Qiagen, Stockach, 
Germany), using the optional on-column DNase treatment. 
Upper and lower jaws were dissected from multiple indi- 
viduals and pooled for each stage, while biological repli- 
cates were kept separate. RNA quality was confirmed 
using agarose gel electrophoresis and purity and quantity 
were assessed with a NanoVue spectrophotometer (GE 
Healthcare, Freiburg, Germany). Prior to cDNA synthesis, 
we confirmed that gDNA contamination was negligible 
(when our samples were amplified alongside positive con- 
trols of comparable concentration, products were detected 
> ten cycles after the positive control). cDNA was synthe- 
sized from up to 200 ng of RNA with Superscript III (Invi- 
trogen, Darmstadt, Germany). Primers were designed to 
our candidate genes using AmplifX 1.5.4 [42], with the re- 
verse primer always in the 3 ' UTR, maximizing the chance 
that they would amplify the intended transcripts, rather 
than orthologues or paralogues (Additional file 2: Table SI). 
Primer concentrations were individually optimized to 
minimize dimers (for all our primer combinations either 
no, or negligible dimer levels were detected) and standard 
curves were generated. 

We performed qRT-PCR reactions with a BioRad ClOOO 
Thermal Cycler, using iQ SYBR green Supermix (Bio-Rad, 
Munich, Germany) according to the manufacturer s in- 
structions. Average CT (avgCT) values were calculated 
across three technical replicates and corrected for PCR- 
efficiency (E = PCR-efficiency; E = 10"^^^'^^^^^; relative 
expression = 1/E^^^^^) of the corresponding primer pair as 
described in [43] and [44]. These raw values of relative 
expression of single individuals were then normalized to 
total RNA input determined with RiboGreen RNA assays 
(Invitrogen, Darmstadt, Germany) according to the manu- 
facturers' instructions. After confirming that the data had a 
normal distribution using the Shapiro-Wilks test, we used 
two-way ANOVA to determine whether jaw tissue or 
developmental stage, or an interaction between the two, 
influenced relative expression of our candidate genes. 

qRT-PCR of candidate genes in the jaws of medaka 

RNA extraction was performed as described for halfbeak. 
In spite of using approximately ten specimens per tissue 
per timepoint, the yield was too low for reliable qRT-PCR 
(cDNA < 1 ng/reaction), thus our samples were amplified 
using the QuantiTect Whole Transcriptome Kit (Qiagen, 
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Stockach, Germany). Initial trials in which two samples 
were amplified using calm2 from the same template (n = 4) 
indicated that relative quantitation (RQ) could not be 
reliably scaled against RNA input quantity (average SD = 
0305). Instead, we developed a set of housekeeping genes 
(HKG), testing previous HKGs used for qRT-PCR on 
cichlid pharyngeal jaws [45], plus ^-actin and gapdh. 
After optimisation with q-base, we selected three HKGs, 
including Ubiquitin fusion degradation l-like, ^-actin 
and gapdh. Normalising RQ against these HKGs (RQ 
candidate/mean RQ HKGs) resulted in a significant re- 
duction in the SD of samples run from the same templates 
(average. SD = 0.141). All templates were then normal- 
ised to these three HKGs, for each primer combination 
(Additional file 2: Table SI). We analysed the expression 
of all genes included in the halfbeak analysis with the 
exception of runx2, which could not be reliably amplified 
in medaka. 

Results 

Characterizing growth heterochronies in halfbeak and 
needlefish 

The relationship between jaw shape and age was examined 
in D. pusilla, O. latipes and B, helone, through plotting 
ratios between lower and upper jaws against scaled age 
using a method modified from [23]. At onset age, all 
three species have equally long upper and lower jaws 
(that is, jaw ratios of close to 1.0) (Figure IB). The 
upper and lower jaws of medaka continue to display 
isometric growth, representing the ancestral condition 
(Figure IB, C). Halfbeak and needlefish both undergo 
accelerated lower jaw growth, as indicated by lower jaw 
to upper jaw ratios (lj:uj) that peak at approximately 1.8 
and 3.0 respectively (Figure IB, C). In needlefish, on the 
other hand, jaw development displays a second phase in 
which uj growth accelerates in relation to Ij until they at- 
tain equal length, a proportion that is then maintained 
throughout their lives (Figure IB, C). Halfbeak lj:uj is then 
maintained at approximately 1.8 through to adulthood 
(data not shown). As needlefish jaw growth includes two 
age-dependent slopes prior to offset age (which is also 
when halfbeak jaws reach their adult ratio), we consider 
needlefish jaw development to represent an example of 
mosaic heterochrony [46]. Ontogenetic scaling analysis 
comparing uj and Ij length demonstrates that the medaka 
Ij grows isometrically with the uj, while halfbeak and 
needlefish display peramorphosis in Ij growth, due to a 
positive shift in the slope of their trajectories in com- 
parison to the ancestral medaka slope (Figure ID, E). 
Additionally, we demonstrate that needlefish jaw growth 
displays an inflection later in development, displaying a 
negative slope that brings the uj and Ij back to an iso- 
metric ratio. This is due to a secondary acceleration of 
uj growth, rather than a reduction in the rate of Ij 



growth, as indicated by allometric growth analyses that 
plot uj and Ij length against log body length (Additional 
file 4: Figure S2). 

Examination of skeletal development 

The developmental basis of lower jaw extension was exam- 
ined using histostaining, comparing needlefish and halfbeak 
development to that of medaka, with a particular focus on 
halfbeak and medaka due to their experimental tractability. 
As skeletal development has previously been described for 
these three species [47-49], aspects of skeletal development 
that related directly to growth heterochrony were the main 
focus of this study. Halfbeak and needlefish at birth and 
hatching respectively, displayed a similar jaw morphology 
to age matched medaka (Figure 2 A and E, data not shown). 
One notable exception is that compared to medaka, the 
Meckel's cartilages of halfbeak and needlefish display a shift 
in their posterior limits, corresponding to a posterior fiision 
of the angulo-articular to the dentaries in halfbeak and 
needlefish [47] (conversely, the angulo-articular is nested 
under the dentaries in medaka). 

Shortly after birth/hatching, halfbeak and needlefish dis- 
play two rostrally-oriented outgrowths of the dentaries at 
their rostro-median contact zones (Figure 2 F, Figure 3C-F, 
Additional file 5: Figure S3). The dentary outgrowths con- 
tinue to extend throughout development and do not bare 
any teeth, even in mature adult stages (Figure 2 F-H, data 
not shown). Hence, they will hereon be referred to as 
'toothless extensions'. At early developmental stages, 
the toothless extensions do not stain with alizarin red, 
and only faintly with alcian blue, which is not visible in 
stained histological sections (Additional file 5: Figure S3). 
The toothless extensions become progressively more ossi- 
fied throughout development in a proximo-distal direction, 
(Figure 2 F-H, Additional file 5: Figure S3), however the 
distal tips remain unossified even in adults (Figure 2G-H, 
data not shown) [47,50]. 

The toothless extensions appear to be morphological 
innovations that are critical to establishing growth het- 
erochrony amongst the Belonoidei. No such structure 
is observed in medaka, whUe they are present in both 
halfbeak and needlefish, and are therefore likely to be 
present throughout the Belonoidei (Figure 3A-F). While 
the teeth on the halfbeak dentary articulate with those of 
the premaxiUa (meaning that there is an otherwise iso- 
metric relationship between uj and Ij), it is the toothless 
extensions that project the lower jaw beyond the anterior 
limit of the premaxilla (Figure 4A-C). 

qRT-PCR profiling of candidate gene expression during 
jaw outgrowth 

The molecular basis of heterochronic growth was examined 
in halfbeak as a first step to determining the molecular 
basis of craniofacial diversity within the Belonoidei. The 
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Figure 2 Comparative siceletal development of medaica, Oryzias latipes and iialfbeal<, Dermogenys pusilla. (A-H) lateral orientation, 
(A-D) 0. latipes, (E-H) D. pusilla, stages are matched based on the relative timing of skeletal ossification. The upper and lower jaws of 
medaka maintain roughly similar proportions throughout development (A-D). At birth, the jaws of halfbeak resemble those of medaka 
(E), however they become elongated throughout development, whereby the lower jaw grows proportionally more than the upper jaw 
(F-H). Abbreviations: den = dentary, dpb = days post birth, dph = days post hatch, e = eye, eth = ethmoid, Ipj = lower pharyngeal jaw, 
mc = Meckel's cartilage, mx = maxilla, pmx = premaxilla. 



expression of candidate genes, including bmp2, bmp4, 
runx2, sox9b and three calmodulin orthologues, denoted 
calmly calm2 and calm3 was compared between upper 
and lower jaws at a range of developmental stages that 
appear to be important for jaw outgrowth (Figure 5A, 
Additional file 6: Figure S4, Additional file 3: Table S2). 
We observed significantly differential expression for four 



of the seven genes examined, based on two-way ANOVAs 
(Table 1). The majority of genes in our study were upregu- 
lated in the lower jaw relative to the upper jaw - amongst 
these were sox9b, runx2 and calmly while bmp2 was 
upregulated in the upper jaw (Figure 5B-D, Figure 6D), As 
both sox9b and runx2 were upregulated in the lower jaw, 
we infer that both chondrogenic (cartilage-forming) and 
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Figure 3 Contribution of the toothless extensions to heterochronic jaw growth in representatives of the Belonoidei. Alcian blue/alizarin 
red staining indicates tliat needlefisli {B. belone) and lialfbeal< (D. pusilla) liave tootliless extensions of tlie lower jaw tliat are not present in medal<a 
(0. lotipes). (A, C, E), lateral orientation, (B, D, F) ventral orientation; age matched representatives, including (A, B) 28 dph medaka, 0. latipes, {C, D) 7 
dpb halfbeak, D. pusillo, {E, F) 6 dph needlefish, B. belone. Abbreviations: bl = body length, e = eye, Ij = lower jaw, Ijb = lower jaw basis, te = toothless 
extensions, uj = upper jaw. 



osteogenic (bone-forming) pathways are required for the 
outgrowth of the lower jaw, presumably the result of their 
enhanced expression in Meckel's cartilages and the den- 
taries respectively (Additional file 3: Table S2). Expression 
of bmp2 was significantly higher in the upper jaw, which 
matches neither runx2 nor sox9b (Figure 5B-D), reducing 
the possibility that it is involved in either bone or cartilage 
development for the stages we examined. 

Interestingly, the three calm genes displayed differential 
patterns of expression, where calml was the only calm 
gene whose expression differed significantly between 
the upper and lower jaws based on two-way ANOVAs 
(Figure 6D-F, Table 1). For calm2 and calmS, expression 
remained equal between the upper and lower jaws 
throughout development (Figure 6E and F). Unlike all 
other genes investigated, the relative expression of calml 
progressively increases in the lower jaw as development 
progresses in a pattern that matches jaw outgrowth 
(see accelerated jaw growth for halfbeak in Figure IB; 
Additional file 7: Figure S5). Moreover, medaka does 
not show significant differences in the expression of 
calmly calm2 or calmS in the jaws during development 
(Figure 6A-C, Table 1), providing further evidence that 
calml may control jaw growth in halfbeak. Additionally, 
none of the investigated genes displayed differential 



expression between upper and lower jaws in medaka. 
We consider calml to be a promising candidate for 
heterochronic growth of the lower jaw in halfbeak and by 
extension, of craniofacial diversification in the Belonoidei. 
All genes included in this analysis were under purifying 
selection (dN/dS ratios < 1), indicating that any regulatory 
evolution displayed by our genes of interest is not accom- 
panied by positive selection in their coding regions. 

Discussion 

We characterized the molecular basis of heterochronic 
shifts that underlie striking craniofacial novelties in some 
lineages of beloniform fishes. Previous studies of hetero- 
chrony in this group merely dealt with the overall jaw 
length without taking into account the skeletal elements 
that contributed to jaw length, or the differential expression 
of genes that control their growth [23,27,51]. Our results 
demonstrate that different rates and directions of ac- 
celeration in jaw growth have shaped the evolutionary 
transitions from medaka to halfbeak, and halfbeak to 
needlefish morphologies respectively. Through normalizing 
developmental stage based on skeletal ossification, our 
heterochrony investigations arrive at a different conclu- 
sion than previous studies [2], classifying needlefish jaw 
development as undergoing two phases of acceleration 
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Figure 4 The toothless extensions form the basis for 
heterochronic shifts in jaw growth amongst the Belonoidei. 

Allometric plots displaying ratios between upper and lower jaw 
growth for (A) medaka, (B) halfbeak, and (C) needlefish as well as 
ratios between upper and lower jaw excluding the toothless 
extensions for (B) halfbeak and (C) needlefish. (A) The jaws of 
medaka, which display the ancestral condition, undergo isometric 
growth (represented by the dashed line). (B) Comparison 
between upper and lower jaw growth for halfbeak indicates 
growth acceleration in favour of the lower jaw, while comparison 
between upper and lower jaw excluding the toothless extensions 
indicates isometric growth (represented by the dashed line). (C) A 
similar pattern to halfbeak is observed in needlefish. Note, 
needlefish growth trajectory is based on only the data points 
prior to growth acceleration of the upper jaw. 



(mosaic heterochrony) rather than hypermorphosis. More- 
over, we identify toothless extensions of the dentaries as 
an important morphological innovation that shapes jaw 
heterochrony within the Belonoidei. calml is identified as 
a potential molecular determinant of this craniofacial 
innovation in the Belonoidei. Interestingly, calm has 
been shown to play a role in establishing jaw (and beak) 
morphology in well-known adaptive radiations such as 
Galapagos finches and cichUd fish [17,22]. Our results 
provide evidence that calm may be an evolutionarily 
widely utilized regulator of craniofacial diversity in the 
vertebrates. 

The toothless extensions of the dentaries of the Belonoidei 
represent an important morphological innovation, which 
is likely to have shaped the craniofacial diversity within 
this group, expanding their evolutionary potential and 
their 'evolvability' [52]. Forming a terminal addition to 
the dentaries, the toothless extensions represent a discrete 
and flexible module that displays weak linkage! one of the 
key features of evolvable systems identified by Kirschner 
and Gerhardt [52]. Moreover, the toothless extensions are 
likely to represent an important innovation amongst the 
Belonoidei due to their presence in both halfbeak and 
needlefish. This flexibility in jaw length appears to have 
enabled the Belonoidei to exploit new ecological niches 
over rapid evolutionary timescales. Indeed, the suborder 
Belonoidei includes almost an order of magnitude more 
species than Adrianichthyoidei (240 and 33 species respect- 
ively), which lack this evolutionary innovation. 

The importance of the toothless extensions as an 
evolutionary innovation lies in the highly flexible nature of 
their development and of their ultimate length in the 
adult, attested by the diversity of jaw phenotypes in the 
suborder Belonoidei. For example, there are even remark- 
able differences in lower jaw length among closely related 
halfbeak species [53,54]. In some halfbeak and flying fish 
genera the lower jaw extension is secondarily resorbed after 
its initial outgrowth in the larval stage [53] and lower jaw 
outgrowth has been entirely lost in most flying fish genera, 
as inferred from a phylogenetic reconstruction of their an- 
cestral states [27,51]. It seems plausible that developmental 
pathways controlling cellular proliferation, differentiation 
and apoptosis in the toothless extensions contribute dir- 
ectly to diversity of jaw phenotype in this group. 

In an effort to better understand the molecular determi- 
nants of jaw elongation in halfbeak, expression of a range 
of skeletogenic genes was analysed in several developmental 
stages during jaw outgrowth. We observed dynamic 
patterns of gene expression for many of these, the majority 
of which indicated upregulation in the lower jaw relative 
to the upper jaw - amongst these were sox9b, runx2 
and calml - while bmp2 was upregulated in the upper 
jaw (Table 1). As both sox9b and runx2 were upregulated 
in the lower jaw, we infer that both chondrogenic 
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Figure 5 Expression of candidate genes during jaw outgrowth in Dermogenys pusilla. (A) Cartoon of tissue sampling for tine expression 
analyses. Upper and lower jaws were excised at the level indicated by the orange and blue colouring respectively, from stages indicated. (B-E) 
Expression of candidate genes (B) sox9b, (C) runx2, (D) bmp2 and (E) bmp4. Relative expression values are scaled so that expression in the lower 
jaw in embryos = 1. Standard error is indicated. Two-way ANOVAs indicated statistically significant differences in expression for sox9b {P< 0.001), 
runx2 (P<0.01) and bmp2 (P<0.05) but not bmp4 (P>0.1) (see Table 1). RQ = relative quantitation. 



Table 1 Results of two-way ANOVA on qRT-PCR of 
skeletogenic genes 

Organism Gene Jaw Age Interaction 

(P-value) (P-value) (P-value) 



Significant P-values are highlighted in bold (P<0.05). 



(cartilage-forming) and osteogenic (bone-forming) path- 
ways are required for the outgrowth of the lower jaw. This 
pattern is likely to be the result of their enhanced expres- 
sion in Meckel's cartilages and the dentaries respectively, 
potentially including the toothless extensions. Future 
studies of spatial expression patterns are required to 
determine if this is the case. Elevated relative expression 
of bmp2 in the upper jaw, when both osteogenesis and 
chondrogenesis are enhanced in the lower jaw, suggests 
that bmp2 does not contribute significantly to lower jaw 
outgrowth in halfbeak. Alternatively, expression of bmp2 
may be strong in the teeth (which are at various stages of 
development at birth) rather than the jaws, as bmp2 has a 
conserved role in teleost dentition [55]. The upper jaw 
samples are composed of a relatively higher proportion of 
tooth germs, as the lower jaw also includes toothless 
extensions, potentially diluting the bmp2 mRNA con- 
centration over the entire sample. 

We consider calml to be a strong candidate for elevated 
outgrowth of the lower jaw of D. pusilla, and by extension, 
of craniofacial diversification in the Belonoidei. Expression 
of calml showed a unique pattern compared to its two 
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Figure 6 Expression oi calm paralogues during jaw outgrowth in medaica and halfbealc. (A-F) Expression of col m paralogues in (A-C) medaka 
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significant difference in expression for colml in halfbeak (P< 0.001), but not medaka (P> 0.1). Neither colm2 nor coImB show differential expression for 
halfbeak or medaka {P > 0.1 ) (see Table 1 ). RQ = relative quantitation. 



paralogues, calm2 and calmS, suggesting that calml may 
have undergone regulatory evolution in the Belonoidei, 
resulting in outgrowth of the lower jaw. Moreover me- 
daka, which has an ancestral jaw morphology, displays 
equal expression of calml, calm2 and calmS in the upper 
and lower jaws, strengthening the association between ele- 
vated relative calml expression and lower jaw outgrowth. 
Previous studies suggested that calm promotes craniofa- 
cial novelties in cichlid and finch through different devel- 
opmental mechanisms, which may shed light on its role in 
halfbeak. For cichlids, calm expression was increased in 
the dentaries of deep-jawed species, with a proposed role 
in increasing osteoblast proliferation [22], while in finches 
with robust beaks, calm expression was increased in the 



distal mesenchyme of the frontonasal prominence, which 
appeared to increase the length of the pre-nasal cartilage 
rod [17]. As the dentary develops through intramembra- 
nous rather than endochondral ossification, we predict that 
similar to cichlid, overexpression of calml in the lower jaw 
of the halfbeak directly increases osteoblast proliferation, 
promoting jaw elongation, however in situ hybridization is 
required to determine whether this is the case. Osteo- 
blast proliferation may be promoted by calm, through 
increased activation of Calcium/calmodulin-dependent 
kinase II (CdKII) [56], which controls the expression of 
osteoblast proliferation and differentiation-related genes 
such as Osterix [57] and dlxS [58]. To date, it is unclear 
whether a single calm orthologue underlies morphological 
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novelties amongst halfbeaks, finches and cichlids, as our 
phylogenetic analyses could not assign clear orthology to 
the finch calm paralogue (data not shown), and the cichlid 
calm sequence has not been made publically available. 

In line with other vertebrates, the halfbeak and medaka 
genomes include multiple calm genes that encode proteins 
with identical amino acid sequences (for a review see [59]; 
see also [60-64]). The retention of such apparently redun- 
dant paralogues in the genome may serve a functional 
purpose, allowing for calm expression to be fine-tuned in 
a spatially and temporally appropriate manner, main- 
taining tight regulation of calcium sensitive pathways. 
For example, studies of calm expression in the rat brain 
demonstrate that they differ spatially [65,66], temporally 
[66,67] and in response to cortical lesions [68] and reser- 
pine, an inhibitor of catecholaminergic neurotransmitters 
[69]. Our results suggest that calml has undergone regu- 
latory evolution after the split between Adrianichthyoidei 
and Belonoidei, leading to its increased expression in the 
lower jaw. We hypothesize that the increased expression 
of calml locally increases the activation of calcium sensitive 
pathways, activating osteoblast proliferation, which leads to 
lower jaw elongation [70-72]. 

Our results suggest that calml expression controls the 
heterochronic growth of a highly modular and develop- 
mentally flexible trait, the toothless extension, potentially 
promoting enhanced trophic differentiation amongst the 
Belonoidei. The variable growth rates of the jaws of Belo- 
noid fishes allows us to test this hypothesis - for example, 
as the ratio between the lower and upper jaws is higher in 
B. belone than in D, pusilla, we predict that the ratio of 
relative expression between lower and upper jaws would 
be higher in B, belone than D, pusilla. It would be particu- 
larly useful to test this hypothesis in beloniform genera that 
display marked differences in the relative lengths of their 
upper and lower jaws, such as Nomorhamphus spp. [54]. 

In their landmark paper on evolvability, Kirschner and 
Gerhardt highlighted the Calmodulin pathway as a particu- 
larly evolvable system due to the flexibility in its binding to 
target proteins [52]. Our study demonstrates a further 
feature of Calmodulin that may promote evolvability: a 
reduced level of pleiotropy of this essential pathway, 
achieved through the retention of duplicated gene copies 
that encode identical proteins. This may explain why Cal- 
modulin has been repeatedly utilised in the development 
of diverse craniofacial innovations [17,19,22]. This allows 
for regulatory evolution to act on individual gene copies, 
while the functions of the other copies remain unaffected. 
The presence of multiple calm gene copies in the genome 
extends the modularity of calm expression through re- 
ducing the pleiotropy of this pathway. This feature of 
Calmodulin may help to explain its repeated implication 
in the generation of morphological novelties amongst 
vertebrates [17,22]. 



Conclusions 

Heterochronic shifts in growth can result in dramatic 
alterations to morphology, which sometimes lead to a 
new adaptive peak via a simple developmental switch. 
Our fine-scale study of growth heterochrony in three 
beloniform fish species identified toothless extensions 
of the dentaries as an important source of heterochro- 
nic growth amongst the Belonoidei. This suborder is 
considerably more speciose than its' sister group the 
Adrianichthyoidei, which has retained the ancestral morph- 
ology. Regulatory evolution of calml is likely to contribute 
to heterochronic growth of the lower jaw amongst the 
Belonoidei, facilitating their impressive diversification 
and promoting their evolutionary success. 

Additional files 



Additional file 1: Figure SI. Examples of heterochronic shifts ir^ growth. 

Additional file 2: Table 51. Primers used ir^ this study. 

Additional file 3: Table 52. Skeletal elemer^ts cor^tributing to ler^gth 
measurements. 

Additional file 4: Figure 52. Further investigation of heterochrony and 
allometry. 

Additional file 5: Figure 53. Histological analysis of jaw development 
in D. pusillo. 

Additional file 6: Figure 54. Phylogenetic analysis of genes cloned 
from D. pusilla. 

Additional file 7: Figure 55. Ratio of gene expression in upper and 
lower jaws of D. pusillo. 



Abbreviations 

bp: base pairs; CdKII: Calcium/calmodulin-dependent kinase II; G": threshold cycle; 
dpb: days post birth; dph: days post hatch; eth: ethmoid; HKG: housekeeping 
gene; Ij: lower jaw; mc: Meckel's cartilage; ML: maximum likelihood; mx: maxilla; 
pmx: premaxilla; RQ: relative quantitation; te: toothless extensions; uj: upper jaw. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

AM and HG designed the study. HG and CK conducted the laboratory work. 
HG wrote the first draft and all authors revised the manuscript, and 
approved it for final publication. 

Authors' information 

Helen Gunter is a postdoctoral fellow in the Meyer lab and in the 
Zukunftskolleg of the University of Konstanz, Germany. She investigates the 
environmental and genomic interactions that shape animal phenotypes, 
focusing primarily on teleost craniofacial development. Claudia Koppermann 
conducted her diploma research on halfbeak development at the University 
of Konstanz. Axel Meyer is an evolutionary biologist interested in the origin 
of diversity at different levels of biological organization. 

Acknowledgements 

We wish to thank Snjezana Kuzir for providing skeletal specimens of Belone belone. 
Eckhard Witten for advice and protocols for teleost histology, and Manfred SchartI, 
for providing medaka. The Zukunftskolleg and the DFG (to HG) and the 
University of Konstanz (to HG and AM) provided funding for this research. The 
undergraduate student helpers and staff of the Tierforschungsanlage (TP A) at 
the University of Konstanz assisted with feeding and maintenance of the fish 
used in this study. We thank Frederico Henning, Ralf Schneider, Shaohua Fan, 
Arne Jacobs, Gerrit Begemann, Shigehiro Kuraku, and Nathalie Feinerfor 
assistance and helpful suggestions on this research. 



Gunter et al. EvoDevo 2014, 5:8 
http://www.evodevojournal.conn/content/5/1/8 



Page 12 of 13 



Author details 

^Department of Biology, Lehrstuhl fur Zoologie und Evolutionsbiologie, 
University of Konstanz, Universitatstrasse 10, 78457 Constance, Germany. 
^Zukunftskolleg, University of Konstanz, Universitatstrasse 10, 78457 
Constance, Germany. 

Received: 12 August 2013 Accepted: 12 December 2013 
Published: 5 February 2014 



References 

1 . Carroll SB: Endless Forms Most Beautiful: The New Science of Evo Devo. 
London: W. W. Norton & Company; 2005. 

2. de Beer GR: Ennbryos and Ancestors. London: Oxford University Press; 1940. 

3. Gould SJ: Ontogeny and Phylogeny. Cambridge: Mass: Belknap Press of 
Harvard University Press; 1977. 

4. Raff RA, Wray GA: Heterochrony: developmental mechanisms and 
evolutionary results. J Evol Biol 1989, 2:409-434. 

5. Rice SH: The analysis of ontogenetic trajectories: when a change in size 
or shape is not heterochrony. Proc Natl Acad Sci USA] 997, 94:907-91 2. 

6. Klingenberg CP: Heterochrony and allometry: the analysis of evolutionary 
change in ontogeny. Biol Rev 1998, 73:79-123. 

7. Smith KM: Heterochrony revisited: the evolution of developmental 
sequences. BiolJ Linn Soc 2001, 73:169-186. 

8. Smith KK Time's arrow: heterochrony and the evolution of development. 
IntJDev Biol 2003, 47:613-621. 

9. Alberch P, Gould SJ, Oster GF, Wake DB: Size and shape in ontogeny and 
phylogeny. Paleobiology 1979, 5:296-317. 

1 0. Raff EC, Raff RA: Dissociability, modularity, evolvability. Evol Dev 2000, 
2:235-237. 

11. Hall BK Developmental processes underlying heterochrony as an 
evolutionary mechanism. Can J Zoo 1 1984, 62:1-7. 

12. Yamamoto Y, Jeffery WR: Central role for the lens in cave fish eye 
degeneration. Science 2000, 289:631-633. 

13. Carleton KL, Spady TC, Streelman JT, Kidd MR, McFarland WN, Loew ER: 
Visual sensitivities tuned by heterochronic shifts in opsin gene 
expression. BMC Biol 2008, 6:22. 

14. Albertson RC, Yan YL, Titus TA, Pisano E, Vacchi M, Yelick PC, Detrich HW, 
Postlethwait JH: Molecular pedomorphism underlies craniofacial skeletal 
evolution in Antarctic notothenioid fishes. BMC Evol Biol 2010, 10:4. 

1 5. Kerney R, Gross JB, Hanken J: Early cranial patterning in the direct-developing 
frog Eleutherodactylus coqui revealed through gene expression. Evol Dev 2010, 
12:373-382. 

16. Hall BK Evo-Devo: evolutionary developmental mechanisms. Int J Dev Biol 
2003, 47(7/8):49 1-496. 

1 7. Abzhanov A, Kuo WP, Hartmann C, Grant BR, Grant PR, Tabin CJ: The 
calmodulin pathway and evolution of elongated beak morphology in 
Darwin's finches. Nature 2006, 442:563-567. 

18. Abzhanov A, Protas M, Grant BR, Grant PR, Tabin CJ: Bmp4 and 
morphological variation of beaks in Darwin's finches. Science 2004, 
305:1462-1465. 

19. Grant PR, Grant BR, Abzhanov A: A developing paradigm for the 
development of bird beaks. Biol J Linn Soc 2006, 88:1 7-22. 

20. Albertson RC, Kocher TD: Genetic and developmental basis of cichlid 
trophic diversity. Heredity 2006, 97:21 1-221. 

21. Albertson RC, Streelman JT, Kocher TD, Yelick PC: Integration and evolution 
of the cichlid mandible: the molecular basis of alternate feeding 
strategies. Proc Natl Acad Sci USA 2005, 102:16287-16292. 

22. Parsons KJ, Albertson RC: Roles for Bmp4 and CaMI in shaping the jaw: 
evo-devo and beyond. Annu Rev Genet 2009, 43:369-388. 

23. Boughton DA, Collette BB, McCune AR: Heterochrony in jaw morphology 
of the needlefishes (Teleostei: Belonidae). SystZool 1991, 40:329-354. 

24. Schlesinger Z: Zur phylogenie und ethologie der Scombresociden. 
VerL) Zool-Bot Ges Wien 1909, 59:302-339. 

25. Regan CT: The classification of the teleostean fishes of the order 
Synentognathi. Ann Mag N Hist 191 1, 7:327-335. 

26. Nichols JT, Breder CMJ: An annotated list of the Synentognathi with 
remarks on their development and relationships. Zoologica (N Y) 1928, 
8:423-448. 

27. Lovejoy NR, Iranpour M, Collette BB: Phylogeny and jaw ontogeny of 
beloniform fishes. Integr Comp Biol 2004, 44:366-377. 



28. Levit GS, Hossfeld U, Olsson L: The integration of Darwinism and 
evolutionary morphology: Alexej Nikolajevich Severtzoff (1866 to 1936) 
and the developmental basis of evolutionary change. J Exp Zool Part B 
2004, 3028:343-354. 

29. Severtzoff AN: Uber die beziehung zwischen der ontogenese und der 
phylogenese der tiere. Jena Zeitrsctir Naturwiss 1927, 63:51-180. 

30. Lovejoy NR: Reinterpreting recapitulation: systematics of needlefishes 
and their allies (Teleostei: Beloniformes). Evolution 2000, 54(4):1 349-1 362. 

31. Parenti L: A phylogenetic analysis and taxonomic revision of ricefishes, 
Oryzias and relatives (Beloniformes, Adrianichthyidae). Zool J Linn 
Soc-Lond 2008, 154:494-610. 

32. Hansen TP: The evolution of genetic architecture. Ann Rev Ecol Evol Syst 
2006, 37:123-157. 

33. Near TJ, Dornburg A, Eytan Rl, Keck BP, Smith WL, Khun KL, Moore JA, Price 
SA, Burbrink FT, Friedman M, et al: Phylogeny and tempo of diversification 
in the superradiation of spiny-rayed fishes. Proc Natl Acad Sci USA 2013, 
110:12738-12743. 

34. Kinoshita M, Murata K, Naruse K, Tanaka M: A Laboratory Manual for Medaka 
Biology 1st edition. Wiley-Blackwell: Ames, Iowa, USA; 2009. 

35. Walker MB, Kimmel CB: A two-color acid-free cartilage and bone stain for 
zebrafish larvae. BiotecL) Histoctiem 2007, 82:23-28. 

36. Gould SJ: Of coiled oysters and big brains: how to rescue the terminology 
of heterochrony, now gone astray. Evol Dev 2000, 2(5):241 -248. 

37. Waterhouse AM, Procter JB, Martin DMA, Clamp M, Barton GJ: Jalview 
Version 2-a multiple sequence alignment editor and analysis workbench. 
Bioinformatics 2009, 25:1 1 89-1 191. 

38. Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel 0: New 
algorithms and methods to estimate maximum-likelihood phylogenies: 
assessing the performance of PhyML 3.0. Syst Biol 2010, 59:307-321. 

39. Hurst LD: The Ka/Ks ratio: diagnosing the form of sequence evolution. 
Trends Genet 2002, 1 8(9):486-487. 

40. Yang Z, PAML4: Phylogenetic Analysis by Maximum Likelihood. Mol Biol 
Evol, 24:1586-1591. 

41. Yang Z, Nielsen R: Estimating synonymous and nonsynonymous 
substitution rates under realistic evolutionary models. Mol Biol Evol 2000, 
17(l):32-43. 

42. Amplifx 1.5.4: http://ifrjr.nord.univ-mrs.fr/AmplifX-Home-page. 

43. Cummings ME, Larkins-Ford J, Reilly CR, Wong RY, Ramsey M, Hofmann HA: 
Sexual and social stimuli elicit rapid and contrasting genomic responses. 
Proc Biol Sci 2008, 275:393-402. 

44. Simon P: Q-Gene: processing quantitative real-time RT-PCR data. 
Bioinformatics 2003, 19:1439-1440. 

45. Gunter HM, Fan S, Xiong F, Franchini P, Fruciano C, Meyer A: Shaping 
development through mechanical strain: the transcriptional basis of 
diet-dependent phenotypic plasticity in a cichlid fish. Mol Ecol 2013, 
22:4516-4531. 

46. Mitteroecker P, Gunz P, Weber G, Bookstein F: Regional dissociated 
heterochrony in multivariate analysis. Ann Anat 2004, 186(5):463-470. 

47. Clemen G, Wanninger AC, Greven H: The development of the dentigerous 
bones and teeth in the hemiramphid fish Dermogenys pusillus 
(Atheriniformes, Teleostei). Ann Anat 1997, 179:165-174. 

48. Kuzir S, Kozaric Z, Gjurcevic E, Bazdaric B, Petrinec Z: Osteological 
development of the garfish {Belone belone) larvae. Anat Histol Embryol 
2009, 38:351-354. 

49. Langille RM, Hall BK Development of the head skeleton of the Japanese 
medaka, Oryzias latipes {Je\eoste\). J Morphol 1987, 193:135-158. 

50. Greven H, Wanninger AC, Clemen G: Dentigerous bones and dentition in 
the hemiramphid fish Dermogenys pusillus (Atheriniformes, Teleostei). 

Ann Anat 1997, 179:21-32. 

51 . Lovejoy N: Phylogeny and ontogeny of needlefishes. Integr Comp Biol 

2000, 42:1268-1268. 

52. Krschner M, Gerhart J: Evolvability. Proc Natl Acad Sci USA] 998, 
95:8420-8427. 

53. Collette BB, McGowen GE, Parin NV, Mito S: Beloniformes: development 
and relationships. Am Soc Ichthyologists and Herpetologists 1984, 1:335-354. 

54. Meisner AD: Phylogenetic systematics of the viviparous halfbeak genera 
Dermogenys and Nomorhamphus (Teleostei: Hemiramphidae: 
Zenarchopterinae). Zool J Linn Soc 2001 , 1 33:1 99-283. 

55. Wise SB, Stock DW: Conservation and divergence of Bmp2a, Bmp2b, and 
Bmp4 expression patterns within and between dentitions of teleost 
fishes. Evol Dev 2006, 8:51 1-523. 



Gunter et al. EvoDevo 2014, 5:8 
http://www.evodevojournal.conn/content/5/1/8 



Page 13 of 13 



56. Zayzafoon M: Calcium/calmodulin signaling controls osteoblast growth 
and differentiation. J Cell Biochem 2006, 97(1):56-70. 

57. Choi YH, Choi J-H, Oh J-W, Lee K-Y: Calmodulin-dependent kinase II regulates 
osteoblast differentiation through regulation of osterix. Biochem Bioph Res 
Co 2013, 432:248-255. 

58. Seo JH, Jin Y-H, Jeong HM, Kim Y-J, Jeong HG, Yeo C-Y, Lee K-Y: 
Calmodulin-dependent kinase II regulates Dlx5 during osteoblast 
differentiation. Biochem Bioph Res Co 2009, 384(1 ):1 00-1 04 

59. Toutenhoofd SL, Strehler EE: The calmodulin multigene family as a unique 
case of genetic redundancy: multiple levels of regulation to provide 
spatial and temporal control of calmodulin pools? Cell Calcium 2000, 
28:83-96. 

60. Nojima H: Structural organization of multiple rat calmodulin genes. J Mol 
Biol 1989, 208:269-282. 

61 . Matsuo K, Sato K, ll<eshima H, Shimoda K, Takano T: Four synonymous genes 
encode calmodulin in the teleost fish, medaka {Oryzias latipes): conservation 
of the multigene one-protein principle. Gene 1992, 119:279-281. 

62. Berchtold MW, Egli R, Rhyner JA, Hameister H, Strehler EE: Localization of 
the human bona fide calmodulin genes CALMI, CALI\/12, and CALI\/13 
to chromosome-14q24-Q31, chromosome-2p21.1-P21.3, and 
chromosome-19q13.2-Q13.3. Genomics 1993, 16:461-465. 

63. Friedberg F, Taliaferro LT: Calmodulin genes in zebrafish (revisited). 
Mol Biol Rep 2005, 32:55-60. 

64. Luan J, Liu ZH, Zhang SC, Li HY, Fan CX, Li L: Characterization, evolution 
and expression of the calmodulini genes from the amphioxus 
Branchiostoma belcheri tsingtauense. Acta Bioch Bioph Sin 2007, 
39:255-264 

65. Berry F, Brown IR: Developmental expression of calmodulin mRNA and 
protein in regions of the postnatal rat brain. J Neurosci Res 1995, 

42:613-622. 

66. Ni B, Landry CF, Brown IR: Developmental expression of neuronal 
calmodulin mRNA species in the rat brain analyzed by in situ 
hybridization. J Neurosci Res 1992, 33:559-567. 

67. Weinman J, Delia Gaspera B, Dautigny A, Pham Dinh D, Wang J, Nojima H, 
Weinman S: Developmental regulation of calmodulin gene expression in 
rat brain and skeletal muscle. Cell Reg u 1 1991, 2:819-826. 

68. Landry CF, Ivy GO, Brown IR: Effect of a discrete dorsal forebrain lesion in 
the rat on the expression of neuronal and glial-specific genes: induction 
of calmodulin, NF-L, SCI, and GFAP mRNA. J Neurosci Res 1992, 
32:280-289. 

69. Ni B, Brown IR: Modulation of a neuronal calmodulin mRNA species in 
the rat brain stem by reserpine. Neurochem Res 1993, 18:185-192. 

70. Whitfield JF, Boynton AL, Macmanus JP, Rixon RH, Sikorska M, Tsang B, 
Walker PR: The roles of calcium and cyclic AMP in cell proliferation. 
Ann NY Acad Sci 1980, 339:216-240. 

71. Godwin SL, Soltoff SP: Extracellular calcium and platelet-derived growth 
factor promote receptor-mediated chemotaxis in osteoblasts through 
different signaling pathways. J Biol Chem 1 997, 272:1 1 307-1 1312. 

72. Yamaguchi T, Chattopadhyay N, Kifor 0, Butters RR, Sugimoto T, Brown EM: 
Mouse osteoblastic cell line (MC3T3-E1) expresses extracellular calcium 
(Ca-0(2-i-))-sensing receptor and its agonists stimulate chemotaxis and 
proliferation of MC3T3-E1 cells. J Bone Miner Res 1998, 13:1530-1538. 



doi:1 0.1 186/2041 -91 39-5-8 

Cite this article as: Gunter et al.: Revisiting de Beer's textbook example 
of heterochrony and jaw elongation in fish: calmodulin expression 
reflects heterochronic growth, and underlies morphological innovation 
in the jaws of belonoid fishes. EvoDevo 2014 5:8. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



